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I In superconductors where the coherence length is comparable to the Fermi 

wavelength, the vortex viscosity depends on the velocity of the vortex, leading 

to non linear equations of motion. The trajectories of vortices driven by 

a. c. fields show a variety of behaviors as function of frequency. Finite 

pertubations give rise to very long lived transients. The relevance of these 

' results to experiments in high-Tc superconductors is discussed. 
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I. INTRODUCTION. 



In the copper oxide superconductors, the coherence length, ^, is similar 
to the Fermi wavelength, /c]^^, and to the lattice spacing, a. The separation 
of energy levels within vortex cores, 6e can be comparable, or larger, than 
the width of these levels (ultra clean limit) and the temperature. Using stan- 
dard parameters, we estimate that 6e ~ 40K, in agreement with experimental 
observations |j^] . The quantization of levels should change drastically dissipa- 
tion processes at temperatures lower that 6e. The standard theory of flux flow 
dissipation breaks down, as it is based on the existence of fast relaxation 
processes within the vortex core. Note that the limit T <^ 5e is also different 
from what is usually defined as the ultra clean limit, in which Tlt <C Se, but 
6e<T§. 

Simple arguments show that energy dissipation due to vortex motion 
should be exponentially reduced at low velocities. The vortex viscosity be- 
haves as r]{v) ~ e"""/" Q, where vq is given by vq ~ 5e^/fi. The exis- 
tence of this threshold velocity can be understood very simply: it corresponds 
to the situation for which the inverse of the vortex-impurity collision time, 
TcoU ~ i/v, becomes comparable to 5e/h. Taking ^ lOA, we find ~ 4 
Km/s. 

The presence of such a term in the energy dissipation of moving complex 
quantum objects with internal levels (atoms, molecules) is rather common 
m. The existence of a gap in the excitation spectrum of the system implies 
that energy dissipation is exponentially suppressed at low temperatures or 
velocities. 

The resulting equations of motion for the vortices become highly non lin- 
ear. In the following, we study the dynamics of vortices, in the pinned regime, 
driven by a. c. fields. Our aim is to contribute to the understanding of re- 
cent experiments which probe the response of vortices at low temperatures 

. These experiments show a viscosity systematically lower than the value 
predicted by the Bardeen-Stephen theory. These experiments have been per- 
formed in the range of 10^ to 10^ Hz. The previous analysis gives t~[^ ~ 10^^ 
Hz. The number of collisions per cycle is large enough to justify the theory 
presented in [^. 

A velocity dependent viscosity has also been analyzed in [Q. In this work, 
it arises from the influence of the pinning centers on an (almost) freely moving 
vortex, at relatively high temperatures. The total viscosity is always finite, 
and its value is comparable to the Bardeen-Stephen prediction. We do not 
consider that regime in the following. 

II. EQUATION OF MOTION. 

We consider pinned vortices, with a linear restoring force. There is some 
controversy on the magnitude of the Magnus force. To keep the number of 
parameters to a minimum, we assume that vortices are carried by the applied 
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current, and that Galilean invariance fully determines this contribution to 
the motion j^. Finally, we describe the external probe through the superfluid 
velocity, with a given time dependence. 

Then, the equations of motion, in the plane perpendicular to the vortex 
axis, are: 

viv)vx = —kpX + avy 

r]{v)vy = -kpy - a[vx - Vs{t)] 

and r]{v) = r]Qe~^'^^^ . We use units such that kp = 1, rjQ = 1 and vq = 1. 
The equations (||) are determined by the ratio a/r]Q. The only difference 
with the standard theory [^-^ lies in the assumption of a velocity dependent 
viscosity. In the presence of a periodic supercurrent, Vs{t) = VgCos{iot), we 
need also to specify the dimensionless values of aw and v^. 

We now solve numerically eqs. (||). We consider first the case of a periodic 
driving current. We choose a = 1 and = 0.5. The value of a is proportional 
to Us, the condensate density, and tjq oc riimp, the concentration of impurities 
Q. Hence, a/rjo ~ ng/riimp- Finally, vq is a fraction of the bulk depairing 
velocity, vo/vdp ~ Se/A, so that should not be greater than vq. 

If we could take t] = 0, the solution of eqs. ([l|) depends only on the 
dimensionless ratio kp/{auj). 

The dynamics of the vortex depend strongly on the value of the driving 
frequency. To a first approximation, we can neglect the friction term in eqs. 
(|l]). Then, the motion of the vortex is purely inductive, and the Hall angle 
changes from 0(ix'<Cl)to7r/2(tt;»l). This picture is qualitatively correct, 
as seen in fig.(^). The unit of length is rj^vo/kp. Taking typical experimental 
values (at low temperatures) for fep ~ 2 x 10^ N / m^ and % ~ 10^^ N s / 
m P,M| , our unit of length is ~ 40A. 

The vortex orbit, however, contains many higher harmonics (see next 
section). At sufficiently high frequencies, these higher harmonics dominate 
the trajectory, and the role of the dissipation needs to be taken into account. 

The corresponding velocities are plotted in fig.(^). In order to interpret 
this result, we need to consider the influence of the dissipation. Firstly, it 
limits the maximum values of the vortex velocity, which never becomes much 
greater than 1, even close to resonance. Above the resonant frequency, the 
motion becomes highly irregular. Note that the number of degrees of freedom 
in eqs.(2) suffices to generate chaotic behavior in the presence of a periodic 
driving force. We have checked that eqs.(||) do, indeed, show chaotic behavior, 
for a <C 7?o and high frequencies. 

We have also analyzed the response to a finite pulse, as in the experiments 
reported in We find very long lived transients after the pulse has been 
switched off. This effect can be understood analytically. In the absence of a 
driving force, we can take the time derivative of eqs. (Q), to obtain: 

ri{v) ^1 + {VxVx + VyVy) = -kpV^ + a{VyVx - VxVy) 

r]{v){VxVy - VyVx) = a{VxVx + VyVy) 



which gives: 
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At long times, v <C and ri{v) <^ a. Then: 



kpvriiv) 



V^e-^ (4) 



a 

which, to logarithmic accuracy, yields: 

hm V ~ ^ — TY (5) 

Our numerical results are consistent with this asymptotic behavior. Note 
that, if the viscosity was independent of velocity, transients decay exponen- 
tially, with a well defined relaxation time, cx ?//(?7^ + a^). 



III. NON LINEAR RESISTIVITY. 

From the vortex velocities, plotted in fig. (^), we infer the induced voltage, 
by means of the Josephson relation, V ^ vxB. Hence, pxx oc Vy and pxy oc Vx- 
In order to analyze numerically the generation of higher harmonics, solve eqs. 
@) in the presence of a periodic current, and perform a Fourier transform. 

The linear resistivities are plotted in fig. The two components which 
are finite in the limit of vanishing dissipation, Im/?xx and Repxy are approxi- 
mately well described by the standard solution of eqs. (0). We find a sharp 
resonance at w ~ kp/a. The resonance, however, is very asymmetric, with a 
sharp rise at low frequencies, and a slow decay at high frequencies. 

The other components of the resistivity, Repxx and Im/Jxy display more 
unusual behavior. Both are very small at low frequencies. They become 
comparable to the other components near the resonance frequency, and slowly 
decrease at high frequencies. Note that, if the effective viscosity remains finite 
as w — > oo, these functions should tend to a constant. At low frequencies, 
we find the scaling law pxy{io) ~ Pxxi^) l@l '^^ function of frequency. The 
eff'ective viscosity near the resonance can be understood from the standard, 
linear, solution to eqs. @. Near the resonance, the vortex velocity behaves 
as f ~ aVg/r]{v). As v cannot exceed 1, we find that the value of r]{v) at 
resonance should be similar to the value of a. 

If the results shown in fig.@ were analyzed in terms of an effective, 
frequency dependent viscosity, we conclude that r]eff ~ for w <C kp/a, 
rjeff ~ min(a,77o) at resonance, uj ~ kp/a, and that ^ at high frequen- 
cies. 

The resistivities also show a significant number of higher harmonics. At 
low frequencies, this effect is more pronounced in Repxx and Impxy; where 
the standard analysis indicates that they vanish as 77 — > 0. It is unclear to 
us the contribution to this effect of the long lived transients discussed in the 
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preceding section. The value of these harmonics is shown in fig (note 
that the scales are different from those in fig (|3|). The harmonics show also 
a maximum near the resonance frequency. Note that the third harmonic is 
greater than the second. In general, we find the odd harmonics to be larger 
that the even ones. 

IV. CONCLUSIONS. 

We have studied the motion of pinned vortices in the limit where the 
spacing between quasiparticle states within the core is greater than the inverse 
scattering time and the temperature. In this regime, the conventional analysis 
of flux flow dissipation ceases to be valid. The dependence of the energy 
dissipation on the vortex velocity is highly non linear. There is a crossover 
velocity, above which dissipation can be regarded as "conventional ", and 
below which dissipation is exponentially suppressed. 

The non linear effects dominate the behavior of Repxx and Impxy at low 
frequencies. There is an asymmetric resonance around oJres ~ kp/a, where 
kp is the pinning constant, a = husjl is the magnitude of the Magnus force 
( Ug is the condensate density). In terms of a frequency dependent viscosity, 
we find that r\{ijj) ~ below the resonance, r](ujres) ~ o), where % is 

comparable to the Bardeen-Stephen value, and that ry(a;) — > as a; — > oo. 

We have not attempted to average over a distribution of pinning constants. 
To a first approximation, it is equivalent to an average over frequencies, in 
a plot like the one shown in fig. (^). The most noticeable effect will be 
the smoothing of the resonance shown there. We have also not analyzed the 
interaction between vortices pinned in different regimes. The voltages induced 
by the vortex motion tend to be higher for those vortices close to to resonance, 
and it is likely that they will entrain the others. Finally, we do not consider 
here additional sources of dissipation, like the existence of subgap states in a 
d-wave superconductor Q. These effects will be analyzed elsewhere. 

The experiments reported so far focus mainly on the temperature depen- 
dence of the vortex viscosity and pinning constant Pulse experiments, 
like those presented in 0, may be difficult to interpret, as the exponential 
suppression of the vortex viscosity at low velocities implies the existence of 
long lived transients (which may be observed, however, as echo signals). It 
would be interesting if the non linear effects discussed here can be confirmed 
experimentally. 
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FIGURES 



FIG. 1. Trajectories followed by vortices driven at different frequencies, a) a; = 0.05. b) 
uj = 0.15. c) (J = 0.5. d) (J = 5. The "ideal " resonance frequency described in the text is, in 
these units, ojres = kp/{2TTa) = 0.159. For typical experimental parameters, u^es ^ 2 x 10^^ Hz. 
As mentioned in the text, a reasonable unit of length is ~ 

FIG. 2. Time dependence of the vortex velocities for the four cases shown in fig. 1. The unit 
of velocity is w o ~ 4 Km/s. Full line, Vy. Dashed line Vx- T is the period. 

FIG. 3. Frequency dependence of the vortex resistivities as function of frequency. The lowest 
harmonic of the resistivity is shown here. Full line. Re p. Dotted line, Im p. Resistivity units are 
arbitrary. 

FIG. 4. Second (upper panels) and third (lower panels) harmonics of the resistivity. The 
conventions are as in fig.3. The units are the same as the ones used in fig. 3. 
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